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Abstract 

Comparative statistical properties of Parkinson, Garman-Klass, 
Roger-Satchell and bridge oscillation estimators are discussed. Point 
and interval estimations, related with mentioned estimators are con- 
sidered 

1 Examples of volatility estimators 

Consider dependence on time t of the price P{t) of some financial instru- 
ment. As a rule, at discussing of volatility, one consider its logarithm 

X{t) := \nP{t). 

Let point out one of the conventional volatility V{T) definition, which we 
are using in this paper: It is the variance 

V{T) := Var [Y{t, T)] = E [Y^t, T)] - [Y{t, T)] . (1) 

of the log-price increment Y{t,T) := X{t + T) — X{t) within given time 
interval duration T. 

Recall, Garman-Klass (G&K) |Tj, Parkinson (PARK) [2] and Roger-Satch- 
ell (R&S) [3] volatility estimators are resting on the high and low values: 

H := sup Y(t,t'), L = inf Y(t,t'). (2) 

t'e{o,T) t'mT) 
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Accordingly, PARK estimator is equal to 

V - - ^'^ (3) 
In 16 • 

while G&K estimator given by expression 

Vg ■= h{H - Lf - k2{C{H - L) - 2HL) - hC\ , . 

A;i = 0.511, A;2 = 0.0109, fcg = 0.383. ^' 

Here C := Y{t,T) is the close value of the log-price increment. Recall else 
R&S estimator, equal to 

Vr:= H{H -C) + L{L-C). (5) 

Besides of mentioned well-known estimators, we discuss bridge oscillation 
estimator. Below we call it shortly by bridge estimator. Before to define it, 
recall bridge Z{t, t') stochastic process definition. It is equal to 

Z{t,t') := Y{t,t') - ^ Yit,T), t' E (0,r). (6) 

Let introduce high and low of the bridge: 

n:= max Z(t,t'), C := min Z(t,t'). (7) 

t'e(o,T) t'e{o,T) 

Accordingly, mentioned above bridge volatility estimator given by 

H ■.= Kin- cf . (8) 

The value of the factor k, will be calculated later. 



2 Geometric Brownian motion 

One of conventional models of price stochastic behavior is geometric 
Brownian motion (see [1H6]). In particular, it is used in theoretical justi- 
fication of G&K, PARK and R&S estimators. Below we discuss statistics 
of mentioned volatility estimators in frame of geometric Brownian motion 
model. Namely, we assume that increment of the log-price is of the form 

Y{t,T) = fiT + aB{T). (9) 

Here /i is the drift of the price, while B{t) is the standard Brownian motion 
B{t) ~ J\f{0,t). Factor o"^ is the intensity of the Brownian motion. 
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Recall, Brownian motion posses by self-similar property 

Bit) Vt B (^^y VT>0, (10) 

where and below sign ~ means identity in law. 

Using pointed out self-similar property, one can ensure that 

Y{t,t') ~ (tVt x(r,7), 

x{t, 7) := 7-r + B{t), 7 := ^VT, r := - G (0, 1). 

a 1 

Henceforth we call process a;(r, 7) by canonical Brownian motion, while factor 
7 by canonical drift. Using relations ([3j), (jl]), ([8]) and (ITT]) , one find that 

~ V{T) ■ {),(7), Vy ~ V{T) ■ {;,(7), H ~ V{T) ■ v^, 

Vrr.V{T)-Vr{i), V{T)=a^T. 

We have used above canonical estimators: 
(f 

^^^^^■^hUe' ^b-=i^s', d:=h-l, s:=C-(, ^^2) 
^9(7) '■— ^1*^^ ~ A;2(c(i — 2hc) — k^c^, Vr = h{h — c) + l{l — c), 
containing high, low and close values 

h := sup x(r,7), / := inf a:(r,7), c:=2;(l,7), (13) 

tG(0,1) ^6(0,1) 

of canonical Brownian motion, and high and low values 

^= sup z{r), C-= inf ^2:(r), (14) 
Te(o,i) re(o,i) 

of the canonical bridge 

z{t) := x(r, 7) - Tx{l, 7) = B{t) - r ■ B{1), r G (0, 1). (15) 

Plots of the typical paths of the canonical Brownian motion x(r, 7) ffTTj) for 
7 = 1 and corresponding canonical bridge z{t) f|T5|) are given in figure [TJ 

It is worthwhile to note that the closer expected values of canonical esti- 
mators ^^(7), "1)3(7), Vr and vi, to unity, the less biased corresponding original 
volatility estimators. Analogously, the smaller variances of canonical estima- 
tors the more efficient original volatility estimators Vp, Vg, Vr and Vb. 

Notice additionally that canonical drift 7 of the canonical Brownian mo- 
tion x{t, 7) (ITT!) IS, cLS cL rule, unknown. Nevertheless, to get some idea about 
dependence on drift /i of bias and efficiency of volatility estimators, we will 
discuss below in details dependence of canonical estimators statistical prop- 
erties on possible values of the factor 7. 
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Figure 1: Typical paths of canonical Brownian motion x{t, 7) flTTj) for 7 = 1 
and corresponding canonical bridge z{r) f|T5|) 

3 Comparative efficiency of PARK and bridge 
estimators 

Resting on, given at Appendix, analytical formulas for probability density 
functions (pdfs) of random variables (IT^ and dH]), we explore in this section 
some atatistical properties of canonical PARK estimator Vp{'~f) and bridge 
one Vb (fT2l). 

Let check, first of all, unbiasedness of canonical PARK estimator. To 
make it, let calculate, with help of pdf qx{S) flA.7p . mean square of oscillation 
d = h — l of the canonical Brownian motion x{t, 7) at the zero canonical drift 
(7 = 0). After simple calculations obtain 



From here and from expression (|T2|) of canonical PARK estimator "0^(7) one 
can see that the following expression is true 



Let find now the factor n at expressions ([8]) and f lT2|) . To make it, calculate 
first of all the mean square of the bridge oscillation. Due to expression (IA.9|) 




(16) 



Efe(7 = 0)] = 1. 
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Figure 2: Plot of canonical PARK estimator Op (7) mean value, as function of 
canonical drift 7. It is seen that with growth of 7 PARK estimator becomes 
more and more biased. Straight line is the plot of canonical bridge Vb, mean 
value 



for the bridge oscillation s fll2l) pdf, one have 



1 TT^ 



m? 6 

m=l 



E[{)fe] = l K=:^ ^b=-^. (17) 



Accordingly, unbiased canonical bridge estimator has the form 

1 . 

The great advantage of the bridge estimator is its unbiasedness for any 
drift. This remarkable property of the pointed out estimator is the conse- 
quence of the fact that bridge Z(t, t') ([6]) and its canonical counterpart z{t) 
don't depend on the drift ^ (canonical drift 7) at all. On the contrary, PARK 
estimator becomes essentially biased at nonzero drift. In figure |2] depicted 
dependence on 7 of canonical PARK estimator expected value, illustrating 
bias of PARK estimator at nonzero drift. Corresponding curve obtained with 
help of analytical expression (]A.6|) for canonical bridge oscillation pdf. 



Let calculate variances of canonical PARK and bridge estimators. After 
substitution into the rhs of expression 



1 r°° 

E[t)p'(7 = 0)] := 5%{5)d5 
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the sum flA.7p for the canonical Brownian motion oscillation pdf qx{S), and 
after summation obtain for 7 = 0: 

nvlil = 0)] = ^ =^ 1-40733. 
Accordingly, variance of canonical PARK estimator Vp is 

VsltIvJO)] = - 1 ~ 0.407. (18) 

As the next step, we calculate variance of canonical bridge estimator Vb 
(fT7|l . Sought variance is equal to 

Var[t;,] := - E[s'] - 1. 
After substitution here, following from (IA.9p . relation 



7^ 

30' 



obtain 

Var[{)b] = ^ - 1 = 0.2. (19) 



Comparing equalities (fT8|) and ( 1T9|) . one can see that variance of bridge esti- 
mator approximately twice smaller than variance of PARK estimator. 

Recall, variance of bridge estimator does not depend on drift. On the 
contrary, variance of PARK estimator essentially depends on the drift. One 
can see it in figure [3l where depicted plot of dependence, on canonical drift 
7, of canonical PARK estimator variance. 

Notice else that bias of some estimator is insignificant only if it is much 
smaller than rms of corresponding estimator, i.e. is small the relative bias: 

,^ E|f.(7)l - 1 p 

/Var[t)(7)] 



Plot of canonical PARK estimator relative bias, as function of canonical drift 
7 depicted in figure HI 

4 Interval estimations on the basis of PARK 
and bridge estimators 



Given at Appendix analytical expressions (IA.6p . flA.7p and (IA.9P for 



canonical Brownian motion and canonical bridge random oscillations pdfs 
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Figure 3: Plots of dependence on 7 of canonical PARK estimator variance. 
Straight line is the variance of canonical bridge estimator 




Figure 4: Plot of relative bias ( 120|) of canonical PARK estimator as function 
of canonical drift 7 



allow us to explore in details probabilistic properties of PARK and bridge 
canonical estimators. Let find, at first, pdfs of mentioned canonical esti- 
mators random values. It is well known from Probabilistic Theory that pdf 
Wp{x; 7) of canonical PARK estimator is expressed through pdf qx{S; 7) (1A.6|) 
of canonical Brownian motion oscillation by the relation 

l^p(x; 7) = (a/ox; 7) , a = In 16. (21) 

Similarly, pdf of canonical bridge estimator is equal to 



W^{x) = s^—q^[^), a = — . (22) 

Here qh{S) (1A.9P is the pdf of canonical bridge oscillation. Plots of canonical 
PARK estimator pdf, for 7 = 0, and pdf of canonical bridge estimator are 
depicted in figure [51 In figure [6] are comparing pdfs of canonical PARK 
estimator, for 7 = 1, and pdf of canonical bridge estimator. It is seen in 
both figures that pdf of canonical bridge estimator is better concentrated 
around its expected value E[{)f,] = 1 than canonical PARK estimator pdf. 

Knowing estimators pdfs, one can produce interval estimations of pos- 
sible volatility values. Consider typical interval estimation: Let V is some 
volatility estimator, equal to 

V = V{T) ■ V. (23) 

Here v is corresponding canonical estimator, while V{T) is the measured 
volatility. One needs to find probability 



F{N) := Pr|l^(T) < N ■ 



that unknown (random) volatility V{T) is not more than times exceeds 
known (measured) volatility estimated value V. It follows from (l23i) that 
following inequalities are equivalent: 

V{T)<N-V ^ v>l/N. 

Last means in turn that sought probability F{N) is expressed through pdf 
of canonical estimator v by the following way: 



POO 

F{N) = Pr{v> l/N} = / W{x)dx. (24) 

Jl/N 

Here W{x) is the pdf of canonical estimator v. 
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Figure 5: Plots of canonical PARK and bridge estimators pdfs, clearly 
demonstrating "probabilistic preference" of bridge estimator in compare with 
PARK one 




Figure 6: Plots of PARK and bridge canonical estimators pdfs for 7 = 1 
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Figure 7: Plots of probabilities Fp{N) and Fb{N) that true volatility is less 
than times exceeds values of PARK and bridge estimators 

Calculations, resting on relations (12T]) . (l22i) . IHM give probability Fb{2) ~ 
0.918 that true volatility is less than twice of given bridge volatility estimator 
value Vfj. It is substantially larger than analogous probability in the case of 
PARK estimator: Fp{2, 7 = 0) ~ 0.813. 

Plots of probabilities F{N) fl24p dependence on the level A^, for PARK 
estimator (in the case of zero drift fi = 0) and for bridge volatility estimator 
are given in figure [71 

5 Comparative statistics of canonical estima- 
tors 

Above, we explored in detail statistical properties of two, PARK and 
bridge estimators. Here we compare their statistics and statistics of another 
well-known volatility estimators: G&K and R&S one. Despite to previous 
chapters, where we have used known analytical expressions for pdfs of canon- 
ical PARK and the bridge estimators, below we use predominantly results of 
numerical simulations. 

Namely, we produce M ^ 1 numerical simulations of random sequences 

n 1 ^ 

a;n(7) :=7;^ + ^5Ze„, n = 0, 1, . . . , A^, 0:0(7) = 0, (25) 

"^^^ n=l 
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where {e^} are iid Gaussian variables ~ A/'(0,1). Notice that stochastic 
process Xni'j) of discrete argument n rather accurately approximates, for 
large ^ 1, paths of canonical Brownian motion x{t, 7) (|TT|) . 

Knowing M iid sequences {xn{'y)} one can find corresponding iid samples 
of pointed out above canonical estimators. Everywhere below we take number 
of iid samples M and discretization number equal to 

A^ = 5■10^ M = 5-10^ 

Plots in figure [S] demonstrate rather convincingly accuracy of numerical sim- 
ulations. In figure M are given two hundred samples of canonical G&K and 
bridge estimators, ensuring "by naked eye" that canonical bridge estimator 
is more efficient than G&K one. 

In figure [10] are given, obtained by numerical simulations, plots of canon- 
ical G&K, PARK, R&S and bridge estimators mean values, illustrating bias 
of G&K and PARK estimators for nonzero canonical drift 7 7^ 0, and actual 
absence of bias for bridge and R&S estimators. 

Eventually, in figure [12] are given plots of probabilities that true volatility 
V{T) is larger than half of corresponding estimator value and less than twice 
of it: 

Pa := Pr 1^ < ^(^) < 2V^| = ^ W{x)dx. (26) 

It is seen that for any 7 mentioned probability is essentially larger for bridge 
estimator, than for G&K, R&S and PARK estimators. 
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Figure 8: Upper panel: Histogram of M samples of canonical bridge esti- 
mator Ofe. Solid line is the plot of canonical bridge estimator's pdf, given by 
analytical expression fl22p . (IA.9p . Dashed line is the pdf of canonical PARK 
estimator for 7 = 0. Lower panel: Histogram of M samples of canonical 
G&K estimator Vg for 7 = 0. Solid line is the plot of the canonical bridge 
estimator pdf. Dashed line is the canonical PARK estimator pdf for 7 = 
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Figure 9: Plots of two hundreds samples of canonical estimators. Up to down 
are samples of G&K, R&S, bridge and PARK estimators. It is seen even by 
"naked eye" that bridge estimator estimates volatility more accurately than 
another mentioned estimators 
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Figure 10: Mean values v of canonical PARK (■), G&K (♦), R&S (*) and 

bridge (A) estimators. Solid lines are theoretical expectations, borrowing 
from figure [2] 




Figure 11: Estimations D of variance of PARK (■), R&S (*), G&K (♦) and 

bridge (A) canonical estimators. Solid lines are plots of theoretical variances, 
borroved from the figure Ei It is seen that for any 7 bridge estimator's 
variance significantly smaller than variances of another mentioned estimators 
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Figure 12: Estimations of probability Pa (ES]) at different 7 values, for PARK 
(■), R&S (-^), G&K (♦) and bridge (A) estimators. Solid lines are results 
of theoretical calculations, resting on formula fl2Bl) 
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A Probabilistic properties of high, low and 
close values 

Here are given pdfs of random variables {h, I, c) (fT3|) and variables (^, () 
(|T4|) . which one need for canonical estimators (|T2|) statistical analysis. Let 
begin with random variable c = x{l, 7). Obviously, its pdf is 



It is easy to show, additionally, that joint pdf qxiVi 7) of high value h f lT^ 
of canonical Brownian motion x(r, 7) and the close value c = a;(l, 7) is equal 
to 






(A.l) 



1] > 0. 
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In turn, pdf of high value h f[T^ 



given by expression 
Qxiv, 7) = 



2 

— exp 



— 'je^^^ erfc 



V + 1 



r]>0. (A.2) 



Let write here exphcit expression for joint pdf qxi^lAiXil) of random 
variables (/i, /, c) ( IT3|) . Using formulas, given at the monograph [7] and in the 
article [8], one might show that pointed out joint pdf given by: 

qx{vJ,x; 7) = fix; 7)^{vJ\x), 

(A.3) 

xeii.T]), h> i < xi(-x). 

Here l(x) is the unit step function, equal to unity for x > and zero other- 
wise. Besides, above there is function 

sivJlx) ■■= 

00 

mlmTimir] - £),x) + [l - rn)T{m{r] - £) + £,x)] , (A.4) 

m=—oo 

We need, at exploring statistical properties of canonical G&K estimator, 
in joint pdf qx{S,X',7) of canonical Brownian motion x(t, 7) (ITTi) oscillation 
d = h — I and the close value c = x(l,7). As it follows from (lA.Sp . flA.4p . 
mentioned pdf is equal to 



mx 



gx(<5,x;7) = 4/(x;7) 

m=-oo t \ r\ 

[m{5 - \x\)[{\x\ + 2m(5)2 - 1] - (m + l)(|x| + 2m5)] e-^^^^\^\+^^\ ^^'^^ 



After integration above joint pdf over all x values obtain pdf qx{^'il) of 
oscillation d: 



— exp 

TT 



m=— 00 

,'(5(5 + 4m(5 + 27) 
2 exp ' 



(2m^r-l-m(2 + 7^)) + 

(l + e2^^)(l + m(2 + 72)) ) - 27(a((5, 7, m) - a((5, -7, m)) 

5 > 0. 



(A.6) 
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Here have used auxiliary function 

a{S, 7, m) := e^""^^ [1 + m(3 + 7(5 + 2m6 + 7))] x 



erf 



2m5 + 7 



erf 



d + 2m<5 + 7 
72 



5 > 0. 



In particular case of zero drift (7 = 0), one get from flA.6p following expression 

qM = 

(2m- 1)2^2^ 



(2m — 1) exp 



(A.7) 



6>0. 



All statistical properties of high and low values (jH]) of canonical bridge 
(IT^ are defined by their two- fold joint pdf qb{j]A)i given by relation 



qb 



{r], i) = m [mJ-'(m(77 — i?)) + (1 — ra)J^{m{r] — I) + £)] , 



(Ai 



Following from here pdf qb{5) of canonical bridge oscillation s = ^ — C, given 
by equahty 



Qj,(5) = 85^m2(4m2(52-3)e 



-2m2<52 



(5 > 0. 



(A.9) 



m=l 
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